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Abstract. We study the efficient approximability aLbasic graph and logic problems in the 
literature when instances are specified hierarchically as in pa] or are specified by 1-dimensional finite 
narrow periodic specifications as in ]5q |. We show that, for most of the problems n considered when 
specified using k-level-restricted hierarchical specifications or fc-narrow periodic specifications, the 
following holds: 

(i) Let p be any performance guarantee of a polynomial time approximation algorithm for 
n, when instances are specified using standard specifications. Then Ve > 0, n has a polynomial time 
approximation algorithm with performance guarantee (I + e)p. 

(ii) n has a polynomial time approximation scheme when restricted to planar instances. 
These are the first polynomial time approximation schemes for PSPACE-hard hierarchically 

or periodically specified problems. Since several of the problems considered are PSPACE-hard, 
our results provide the first examples of natural PSPACE-hard optimization problems that have 
polynomial time approximation schemes. This answers an open question in Condon et. al. 
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1. Introduction and motivation. Many practical applications of graph theory 
and combinatorial optimization in CAD systems, mechanical engineering, VLSI design 
and software engineering involve processing large objects constructed in a systematic 
manner from smaller and more manageable components. An important example of 
this occurs in VLSI technology. Currently, VLSI circuits can consist of millions of 
transistors. But such large circuits usually have a highly regular design and con- 
sequently are defined systematically, in terms of smaller circuits. As a result, the 
graphs that abstract the structure and operation of the underlying circuits (designs) 
also have a regular structure and are defined systematically in terms of smaller graphs. 
Methods for describing large but regular objects by small descriptions are referred to 
as succinct specifications. Over the last twenty years several theoretical models have 
been put forward to succinctly represent objects such as graphs and circuits, (see for 
example §, |, [□], [l§ [|| |26j |?], |§ g|, Hj, Here, we study two kinds of succinct 



specifications, namely, hierarchical and periodic specifications. 

Hierarchical specifications allow the overall design of an object to be partitioned 
into the design of a collection of modules; which is a much more manageable task 
than producing a complete design in one step. Such a top down (or hierarchical 
design) approach also facilitates the development of computer aided design (CAD) 
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systems, since low-level objects can be incorporated into libraries and can thus be 
made available as submodules to designers of large scale objects. Other areas where 
hierarchical specifications have found applications are VLSI design and layout JIM [l^, 
|55| , finite element analysis, software engineering and datalog queries (see and 
the references therein). Periodic specifications can also be used to define large scale 
systems with highly regular structures. Using periodic specifications, large objects 
are described as repetitive connections of a basic module. Frequently, the modules 
are connected in a linear fashion, but the basic modules can also be repeated in two 
or higher dimensional patterns. Periodic specifications are also used to model time 
variant problems, where the constraints or demands for any one period is the same 
as those for preceding or succeeding periods. Periodic specifications have applications 
in such diverse areas as transportation planning p8| , [43| , p8[ , parallel programming 
@ |§ and VLSI design @ |§. 

Typically, the kinds of hierarchical and periodic specifications studied in the lit- 
erature are generalizations of standard specifications used to describe objects. An 
important feature of both these kinds of specifications is that they can be much more 
concise in describing objects than standard specifications. In particular, the size of 
an object can be exponential in the size of its periodic or hierarchical specifications. 
As a result of this, problems for hierarchically and periodically specified inputs often 
become PSPACE-hard, NEXPTIME-hard, etc. 

In this paper, we concentrate our attention on 



1. the hierarchical specifications of Lengauer 31, |3J, |35| (referred to as L- 
specifications) and 

2. the 1-dimensional finite periodic specifications of Gale and Wankc |58|| 
(referred to as 1-FPN-specifications). 

Both of these specifications have been used to model problems in areas such as CAD 
systems and VLSI design j|5|, [32|, [36| , transportation planning |l(J , parallel program- 
ming p8|, etc. We give formal definitions of these specifications in and §|^. 

Let LI be a problem posed for instances specified using standard specifications. 
For example, if II is a satisfiability problem for CNF formulas, the standard specifica- 
tion is sets of clauses, with each clause being a set of literals. Similarly if II is a graph 
problem, the adjacency matrix representation or the adjacency list representation of 
the edges in the graph are standard specifications. For the rest of the paper, we use 
1. L-LI to denote the problem II, when instances are specified using the hierar- 



chical specifications of Lengauer |35j (see Definition 4.1), and 



2. 1-FPN-II to denote the problem II, when instances are specified using the 



1-dimensional finite periodic specifications of Wanke |58|] (see Definition |5.l|) . 
Thus for example, L-3SAT denotes the problem 3SAT when instances are specified 
using L-specifications and 1-FPN-3SAT denotes the problem 3SAT when instances are 
specified using 1-FPN-specifications. For the rest of this paper, we use the term 
succinct specifications to mean both L-specifications and 1-FPN-specifications. 

2. Summary of results. In this paper, we discuss a natural syntactic restriction 
on the L-specifications and call the resulting specifications level-restricted specifi- 
cations. (For 1-FPN-specifications our notion of level-restricted specifications closely 
coincides with Orlin's notion of narrow specifications |Q.) Most of the problems 
considered in this paper are PSPACE-hard even for level-restricted specifications (see 
||, H, H). Consequently, we focus our attention on devising polynomial time ap- 
proximation algorithms for level restricted L- or 1-FPN-specified problems. Recall 
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that an approximation algorithm for a minimization problem[J IT provides a perfor- 
mance guarantee of p if for every instance / of II, the solution value returned by the 
approximation algorithm is within a factor p of the optimal value for /. A polyno- 
mial time approximation scheme (PTAS) for problem LI is a family of algorithms 
such that, for e > 0, given an instance / of II, there is a polynomial time algorithm 
in the family that returns a solution which is within a factor (1 + e) of the optimal 
value for i\ The main contributions of this paper include the following. 

(i) We design polynomial time approximation algorithms (for arbitrary in- 
stances ) and approximation schemes (for planar instances) for a variety of natural 
PSPACE-hard problems specified using level-restricted L- or 1-FPN-specifications. 
These are the first polynomial time approximation schemes in the literature for "hard" 
problems specified using either L- or 1-FPN-specifications. To obtain our results we 
devise a new technique called the partial expansion. The technique has two desirable 
features. First, it works for a large class of problems and second, it works well for 
both L-specified and 1-FPN-specified problems. 

(ii) For problems specified using level-restricted L- or 1-FPN-specifications, we 
devise polynomial time approximation algorithms with performance guarantees that 
are asymptotically equal to the best possible performance guarantees for the corre- 
sponding problems specified using standard specifications. 

(iii) The results presented in this paper are a step towards finding sufficient 
syntactic restrictions on the L- or 1-FPN-specifications that allow us to specify a 
number of realistic designs in a succinct manner while making them amenable for 
rapid processing. 

Our results provide the first examples of natural PSPACE-complcte problems 
whose optimization versions have polynomial time approximation schemes. Thus they 
affirmatively answer the question posed by Condon, Feigenbaum, Lund and Shor || 
of whether there exist natural classes of PSPACE-hard optimization problems that 
have polynomial time approximation schemes. 

2.1. The meaning of approximation algorithms for succinctly specified 
problems. When objects are represented using L- or 1-FPN-specifications, there are 
several possible ways of defining what it means to "design a polynomial time approxi- 
mation algorithm" . Corresponding to each decision problem LI, specified using either 
L- or 1-FPN-spccifications, we consider four variants of the corresponding optimiza- 
tion problem. We illustrate this with an example. 

Example 1: Consider the minimum vertex cover problem, where the input is an 
L-specification of a graph G. We provide efficient algorithms for the following versions 
of the problem. 

1. The construction problem: Output an L-specification of the set of vertices 
in the approximate vertex cover C. 

2. The size problem: Compute the size of the approximate vertex cover C for 

G. 

3. The query problem: Given any vertex v of G and the path from the root 
to the node in the hierarchy tree (see §2 for the definition of hierarchy tree) in which 
v occurs, determine whether v belongs to the vertex cover C. 

4. The output problem: Output the approximate vertex cover C. 

Note that our algorithms for the four variants of the problem apply to the same 
vertex cover C. Our algorithms for (1), (2) and (3) above run in time polynomial 



A similar definition can be given for maximization problems. 
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in the size of the L-specification rather than the size of the graph obtained by 
expanding the L-specification. Our algorithm for (4) runs in time linear in the size 
of the expanded graph but uses space which is only polynomial in the size of the 
L-specification. □ 

Analogous variants of approximation algorithms can be defined for problems spec- 
ified using 1-FPN-specifications. Therefore, we omit this discussion. 

These variants are natural extensions of the definition of approximation algo- 
rithms for problems specified using standard specifications. This can be seen as 
follows: When instances are specified using standard specifications, the number of 
vertices is polynomial in the size of the description. Given this, any polynomial time 
algorithm to determine if a vertex v of G is in the approximate minimum vertex cover 
can be easily modified to obtain a polynomial time algorithm that lists all the vertices 
of G in the approximate minimum vertex cover. Thus in the case when inputs are 
specified using standard specifications, (3) can be used to solve (2) and (4) in poly- 
nomial time. The above discussion also shows that given an optimization problem 
specified using standard specifications, variants (1), (3) and (4) discussed above are 
polynomial time inter-reducible. 

The approximation algorithms given in this paper have another desirable feature. 
For an optimization problem or a query problem, our algorithms use space and time 
which is a low level polynomial in the size of the hierarchical or the periodic specifica- 
tion. This implies that for graphs of size N, that are specified using specifications of 
size 0(polylog N), the time and space required to solve problems is only 0(polylog N). 
Moreover when we need to output the subset of vertices, subset of edges, etc. corre- 
sponding to a vertex cover, maximum cut, etc., in the expanded graph, our algorithms 
take essentially the same time but substantially less (often logarithmically less) space 
than algorithms that work directly on the expanded graph. The graphs obtained by 
expanding hierarchical or periodic descriptions are frequently too large to fit into the 
main memory of a computer jnj. This is another reason for designing algorithms 
which exploit the regular structure of the underlying graphs. Indeed, most of the 
standard algorithms in the literature assume that the input completely resides in the 
main memory. As a result, even the most efficient algorithms incur a large number of 
page faults while executing on the graphs obtained by expanding the hierarchical or 
periodic specifications. Hence, algorithms designed for solving problems for graphs or 
circuits represented in a standard fashion are often impractical for succinctly specified 
graphs. We refer the reader to ]3l], for more details on this topic. 

The rest of the paper is organized as follows. Section ^ contains discussion of 
related research. In §|| and §^ we give the basic definitions and preliminaries. 
In §|?] we discuss our approximation algorithms for L-specified problems and 1-FPN- 
specified problems. Finally in we give concluding remarks and directions for future 
research. 

3. Related research. In the past, much work has been done on character- 
izing the complexity of various problems when instances are specified using L- or 
1-FPN-specifications. For periodically specified graphs, several researchers || [?]. |T^. 
p6| , [27], |5(| have given efficient algorithms for solving problems such as deter- 
mining STRONGLY CONNECTED COMPONENTS, TESTING FOR EXISTENCE OF CYCLES, 
FINDING MINIMUM COST PATHS BETWEEN A PAIR OF VERTICES, BIPARTITENESS, 

planarity and minimum COST spanning forests. Orlin @ and Wanke HI dis- 
cuss NP- and PSPACE-hardness results for infinite and finite periodically specified 
graphs. 
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For L-specified graphs, Lengauer et al. |32[ |34|, Hi an d Williams et al. |59|| 
have given efficient algorithms to solve several graph theoretic problems including 
2-COLORING, MINIMUM SPANNING FORESTS and planarity testing. Lengauer and 
Wagner j3?[] show that the following problems are PSPACE-hard when graphs are 
L-specified: 3 COLORING, Hamiltonian circuit and path, monotone circuit 

VALUE PROBLEM, NETWORK FLOW, ALTERNATING GRAPH ACCESSIBILITY and MAX- 
IMUM independent set. In Q, Lengauer and Wanke consider a more general 
hierarchical specification of graphs based on graph grammars and gave efficient al- 
gorithms for several basic graph theoretic problems specified using this specification. 
We refer the reader to (l^, |^| for a detailed survey of the work done in the area of 
hierarchical and periodic specifications. 

A substantial amount of research has been done on finding polynomial time ap- 
proximation algorithms with provable worst case guarantees for NP-hard problems. In 
contrast, until recently little work has been done towards investigating the existence 
of polynomial time approximation algorithms for PSPACE-hard problems. As a step 
in this direction, in jfO| |d| we have investigated the existence and non-existence of 
polynomial time approximations for several PSPACE-hard problems for L-specified 
graphs. In [po[ , we considered geometric intersection graphs defined using the hi- 
erarchical specifications (HIL) of Bentley, Ottmann and Widmayer ||. There, we 
devised efficient polynomial time approximation schemes for a number of problems 
for geometric intersection graphs, specified using a restricted form of HIL. 

Condon, et al. (|, ^| also studied the approximability of several PSPACE-hard op- 
timization problems. They characterize PSPACE in terms of probabilistically check- 
able debate systems and use this characterization to investigate the existence and 
non-existence of polynomial time approximation algorithms for a number of basic 
PSPACE-hard optimization problems. 

4. The L-specifications. This section discusses the L-specifications. The fol- 
lowing two definitions are essentially from Lengauer |3^, |35|, |37j . 

Definition 4.1. An L- specification T = (G\, G n ) of a graph is a sequence of 
labeled undirected simple graphs Gi called cells. The graph Gi has rrii edges and Hi 
vertices, pi of the vertices are called pins. The other (rii—pi) vertices are called inner 
vertices, ri of the inner vertices are called nonterminals. The (rij — r{) vertices are 
called terminals. The remaining ni ~ pi — ri vertices of Gi that are neither pins nor 
nonterminals are called explicit vertices. 

Each pin of Gi has a unique label, its name. The pins are assumed to be numbered 
from 1 to pi. Each nonterminal in Gi has two labels (v,t), a name and a type. The 
type t of a nonterminal in Gi is a symbol from G\, Gj_i. The neighbors of a 
nonterminal vertex must be terminals. If a nonterminal vertex v is of the type Gj in 
Gi, then v has degree Pj and each terminal vertex that is a neighbor of v has a distinct 
label (v,l) such that 1 < I <Pj- We say that the neighbor of v labeled (v,l) matches 
the Ith pin of Gj . 

Note that a terminal vertex may be a neighbor of several nonterminal vertices. 
Given an L-specification r, N = X)i<i<n ni denotes the vertex number, and M = 
J2i<i<n mi denotes the edge number of T. The size of T, denoted by size(T), is 
N + M. 

Definition 4.2. Let T = (Gi, G n ) be an L- specification of a graph E{T) and 
let Ti — (Gi, Gi) . The expanded graph E(Y) (i.e. the graph associated with T) is 
obtained as follows: 
k = 1 : E(T) = Gi. 
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k > 1 : Repeat the following step for each nonterminal v of Gk, say of the type Gj: 
delete v and the edges incident on v. Insert a copy of E(Tj) by identifying the I th 
pin of E(Tj) with the node in Gk that is labeled (v,l). The inserted copy of E(Tj) is 
called a subcell of Gk- 

Observe that the expanded graph can have multiple edges although none of the 
Gi have multiple edges. Here however, we only consider simple graphs, i.e. there is 
at most one edge between a pair of vertices. This means that multi edges are treated 
simply as single edges. We assume that T is not redundant in the sense that for each 
j, 1 < j < n, there is a nonterminal v of type Gi in the definition of Gj, j > i. 

The expansion E(T) is the graph associated with the L-specification T with vertex 
number N. For 1 < i < n, Ti = (Gi, Gi) is the L-specification of the graph E(Ti). 
Note that the total number of nodes in E(T) can be 2™ N \ (For example, a complete 
binary tree with 2 n< - N > nodes can be specified using an L-specification of size O(N).) 
To each L-specification T — {G\, G n ), (n > 1), we associate a labeled rooted 
unoriented tree HT(T) depicting the insertions of the copies of the graphs E(T j) 
(1 < j < n — 1), made during the construction of E(T) as follows: (see Figure 4.1) 

Definition 4.3. Let T = {G\, ...,G n ), (n > 1) be an L-specification of the graph 
E(T). The hierarchy tree of T, denoted by HT(T), is the labeled rooted unordered 
tree defined as follows: 

1. Let r be the root of HT(T). The label of r is G n . The children of r in HT(T) 
are in one-to-one correspondence with the nonterminal vertices of G n as follows: The 
label of the child s of r in HT(T) corresponding to the nonterminal vertex (v,Gj) of 
G n is (v, Gj). 

2. For all other vertices s of HT(T) and letting the label of s = (v, Gj), the chil- 
dren of s in HT(T) are in one-to-one correspondence with the nonterminal vertices of 
Gj as follows: The label of the child t of s in HT(T) corresponding to the nonterminal 
vertex (w,Gi) ofGj is (w,G\). 

Given the above definition, we can naturally associate a hierarchy tree corre- 
sponding to each T i, \ < i < n. We denote this tree by HT(Ti). Note that, each 
vertex v of E(T) is either an explicit vertex of G n or is the copy of some explicit vertex 
v' of Gj (1 < j < n) in exactly one copy CJ of the graph E(Tj) inserted during the 
construction of E(T). This enables us to assign v of E(T) to the unique vertex n v of 
the HT(T) given by 

1. if v is a terminal vertex of G n , then n v is the root of HT(T), and 

2. otherwise, v belongs to the node n v that is the root of the hierarchy tree 
HTiTj), corresponding to CJ . 

Given HT(T), the level number of a node in HT(T) is defined as the length of 
the path from the node to the root of the tree. 

As noted in |3f| , L-specifications have the property that for each copy (instance) 
of a nonterminal, a complete boundary description has to be given. Thus if a non- 
terminal has a lot of pins, copying it is costly. Another property of the definition of 
L-specifications is that nonterminals are adjacent only to terminals. These properties 
ensure that the size of the "frontier" (or the number of neighbors) of any nonter- 
minal is polynomial in the size of the specification. These properties weaken the 
L-specifications with respect to other notions of hierarchy involving a substitution 
mechanism that entails implicit connections to pins at a cell boundary 0, |5^] . As a 
result^regular structures such as grids cannot be specified using small L-specifications. 
(see |32| ) ■ In contrast the graph glueing model of Galperin allows a hierarchical 
description of pins; thus the size of the frontier can be exponentially large. As a result, 



APPROXIMATION ALGORITHMS 



7 



graphs such as grids can be represented using descriptions of logarithmic size. How- 
ever as demonstrated in jy], |34|, |57j , these properties seem to be a prerequisite 
for the construction of efficient exact algorithms for L-specified problems. As subse- 
quent sections show, these restrictions are also necessary in part for devising efficient 
approximation algorithms for L-specified problems. The size of the frontier also has 
a significant impact on the complexity of several basic succinctly specified problems. 
For example, several basic NP-hard problems become PSPACE-hard when specified 
using L-specifications (see |37], Q). In contrast, in a recent paper we show that these 
problems typically become NEXPTIME-hard when specified using the graph glueing 



specifications of [jn) (see p5 



By noting Definition 4.1 



it follows that an L-specification is a restricted form of 
a context-free graph grammar. The substitution mechanism glues the pins of cells to 



neighbors of nonterminals representing these cells, as described in Definition 4.2 . Such 
graph grammars are known as hyperedge replacement systems Jl5| or cellular graph 
grammars pS| ]. Two additional restrictions are imposed on cellular graph grammars 
to obtain L-specified graphs. First, for each nonterminal there is only one cell that 
can be substituted. Thus there are no alternatives for substitution. Second, the 
index of the substituted cell has to be smaller than the index of the cell in which 
the nonterminal occurs. The acyclicity condition together with the "no alternatives" 
condition implies that an L-specification defines a unique finite graph. We observe 
that HT(T) is the parse tree of the unique graph generated by the context-free graph 
grammar T. 



Example 2: Figure |4- 1| depicts th e L- specification G = (Gi,G2,G3) and the asso- 



ciate hierarchy tree HT{G). Figure 4.2 depicts the graph E(G) specified by G. The 
correspondence between pins of Gj and neighbors of Gj in G^, j < i, is clear by the 
positions of the vertices and the pins. □ 

4.1. Level-restricted specifications. We discuss level restricted L-specifications 
now. This is also discussed in 

Definition 4.4. An L-specification T — (Gi, G n ) (n > 1), of a graph G is 
1-level-restricted, if for all edges (u,v)of E(T), either 

(1) n u and n v are the same vertex of HT(T). or 

(2) one of n u or n v is the parent of the other in HT(T). 

Extending the above definition we can define fc-level-restricted specifications. An 
L-specification T = {G\, ...,G n ), (n > I), of a graph E(T) is k -level-restricted, if for 
all edges (u, v) of E(T), either 

(1) n u and n v are the same vertex of HT(T) or 

(2) one of n u or n v is an ancestor of the other in HT(T) and the length of the path 
between n u and n v in HT(T) is no more than k. 

We note that for any fixed k > 1, fc-level-restricted L-specifications can still lead 
to graphs that are exponentially large in the sizes of their specifications. Moreover, 
L-specifications (see [[30, 34 , 3^]) for several practical designs are k-level-restricted for 



small values of k. (For example, it is easy to define a complete binary tree with 
2 f2 ( Ar ) nodes by a 4-level-restricted L-specification of size O(N). Note however, that 
the specification depicted in Figure 4.4 is not 4-level-restricted.) For the rest of the 
paper, given a problem LT specified using standard specifications, we use 4-L-LI to 
denote the problem specified using 4-level-restricted L-specifications and fc-L-LI to 
denote the problem specified using fc-level-restricted L-specifications. 

5. 1-FPN-specifications. Next, we give the definition of 4-dimensional peri- 
odic specifications due to Orlin [j48|, Wanke p8[ and Hofting and Wanke p9||. For 
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Fig. 4.1. An L- specification G of a graph E(G), and the associated hierarchy tree HT(G). 
The mapping between the pins and its neighbors is clear by the relative positions of the pins and its 
neighbors. 



the rest of the paper N and Z denotes the set of non-negative integers and integers 
respectively. 

Definition 5.1. Let G(V,E) (referred to as a static graph) be a finite directed 
graph such that each edge (u, v) has an associated non-negative integral weight t UjV . 
The undirected one way infinite graph G°°(V',E') is defined as follows: 

1. V = {v(p) \veV and p G N} 

2. E' = {(u(p), v{p + t u ,v)) ( u iV) G E , t UfV is the weight associated with the 
edge (u, v) and p G N} 

A 1- dimensional periodic specification T (referred to as 1-P- specification) is given by 
r = (G(V,E)) and specifies the graph G°°(V , E') (referred to as 1-P-specified graph). 

A 1-P -specification T is said to be narrow or ldevel-restricted if\f(u,v) G E, 
t u ,v G {0,1}. This implies that V(u(p),«(g)) G E' , \p — q\ < 1. Similarly, a 1-P- 
specification is fc-narrow or fc-level-restricted if\/(u,v) G E, t UtV G {0, 1, . . .k}. 

We note that if we replace N by Z in Definition ^T], we obtain a two way infinite 
periodically specified graph defined in Orlin p8[ . It is sometimes useful to imagine 
a narrow periodically specified graph G°° as being obtained by placing a copy of the 
vertex set V at each integral point (also referred to as lattice point) on the X-axis (or 
the time line) and joining vertices placed on neighboring lattice points in the manner 
specified by the edges in E. 
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Fig. 4.2. The graph E(G) represented by G specified in Figure 4-1- 





G 



G ' 



Fig. 4.3. A static graph G, and the graph G 4 specified by the 1-FPN '-specification V = (G,4). 



G m is the subgraph of the infinite periodic graph G°° induced by the vertices 
associated with nonnegative lattice points less than or equal to m. Formally, 

Definition 5.2. Let G(V,E) denote a static graph. Let G°°{V',E') denote 
the one way infinite 1-PN-specified graph as in Definition 5.1. Let m > be an 
integer specified using binary numerals. Let G m (V m , E m ) be a subgraph ofG°°{V', E') 
induced by the vertices V m = {v(p)\v € V and < p < to}. A 1- dimensional 
finite periodic specification T (referred to as 1-FPN -specification) is given by V = 
{G{V 1 E),m) and specifies the graph G m (referred to as 1-FPN -specified graph). 

An example of a 1-FPN-specified graph appears in Figure O. In psfl , Orlin 
defined the concept of two way infinite 1-dimensional periodically specified 3CNF 
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formulas and the associated 3SAT problem |12|. It is straightforward to restrict Or- 



lin's definition along the lines of Definition 5.1 to define 1-FPN-specified satisfiability 



problems. As a consequence, we omit the definition here. (See [|44|, |48|, |3J] for formal 
definitions of periodically specified satisfiability problems.) We only give an example 
of 1-FPN-specified 3CNF formula to illustrate the concept. 

Example 3: Let U = {2:1,2:2,2:3} be a set of static variables. Let C be a set of 
static clauses given by (xi(0) + x 2 (0) + 2:3(0)) A (2:1(1) + 2:3(0)) A (2:3(1) + x 2 (0)). 
Let F = (U,C,3) be a 1-FPN-specification. Then F specifies the 3 CNF formula 
F 3 {U 3 ,C 3 ) given by 

(a:i(0) + i 2 (0) + 2:3(0)) A ( Xl (l) + 2:3(0)) A (2:3(1) + 2: 2 (0)) /\ 
(2;i(l) + 2:2(1) + 2:3(1)) A (2;i(2) + 2:3(1)) A (2:3(2) + 212(1)) /\ 

(2;i(2) + 2:2(2) + 2:3(2)) A (2;i(3) + 2:3(2)) A (2:3(3) + 2:2(2)) /\ 

(2;i(3) + 2; 2 (3) + 2:3(3)) 

6. Other preliminaries. Recall that a graph is said to be planar if it can be laid 
out in the plane in such a way that there are no crossovers of edges. For the rest of the 
paper, we use l-pl-II, 1-l-pl-LT and 1-fpn-pl-IT to denote the problem II restricted 
to L-specified planar instances, 1-level-restricted L-specified planar instances and 1- 
FPN-specified planar instances respectively. As shown in Lengauer |35|, given an 
L-specification T, there is a polynomial time algorithm to determine if E(T) is planar. 
Similarly as pointed out in |jl8), given a 1-FPN-specification T, there is a polynomial 
time algorithm to determine if E(T) is planar. Thus for solving L- or 1-FPN-specified 
problems restricted to planar instances, we can assume without loss of generality that 
the inputs to our algorithms consist of planar instances. 

Next, we define the problems MAX SAt(S). The definition is essentially an exten- 
sion of the definition of SAt(S) given in Schaefer fl56| . 

Definition 6.1. (Schaefer j^/ J 
Let S = {Ri, i?2j • • • j Rm} be a finite set of finite arity Boolean relations. (A Boolean 
relation is defined to be any subset of {0, 1} P for some integer p > 1. The integer p 
is called the arity of the relation.) An S-formula is a conjunction of clauses each of 
the form i?i(£i, £2, ■ ■ where £1, £2, • • • are distinct, unnegated variables whose number 
matches the arity of Ri,i G {1, • • • m} and Ri is the relation symbol representing the 
relation Ri. The S- satisfiability problem is the problem of deciding whether a given 
S-formula is satisfiable. 

Given a S-formula F, the problem MAX SAt(S) is to determine the maximum number 
simultaneously satisfiable clauses in F. 

As in Schaefer |36|, given S, Rep(S) is the set of relations that are representable 
by existentially quantified S-formulas with constants. 

Recall from [^] that a S-formula / is said to be planar if its associated bipartite 
graph is planar. The problem PL-3SAT ||{J is the problem of determining if a given 
planar 3CNF formula is satisfiable. Lichtenstein |39] showed that the problem pl- 
3 SAT is NP-complete. 
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Next, we define L-specified S-formulas. Such formulas are built by denning larger 
S-formulas in terms of smaller S-formulas. Just as L-specifications of graphs can repre- 
sent graphs that are exponentially larger than the specification, L-specified S-formulas 
can specify formulas that are exponentially larger than the size of the specification. 

Definition 6.2. An instance F = (F 1 (X 1 ),...,F n - 1 (X n - 1 ),F n (X n )) of L- 
SAt(S) is of the form 

F i (X i ) = ( f\ 
i<j<h 

for 1 < i < n where fi are S-formulas. X n = f>, X % , X -, Z\ Zj, 1 < i < n — 1, 
are vectors of Boolean variables such that X 1 - C X 1 , Zj C Z l , < ij < i. Thus. 
Fi is just a S-formula. An instance o/l-SAt(S) specifies a S-formula E(F), that is 
obtained by expanding the Fj, 2 < j < n, where the set of variables Z's introduced in 
any expansion are considered distinct. The problem L-SAt(S) is to decide whether the 
formula E(F) specified by F is satisfiable. The corresponding optimization problems 
denoted by L-MAX-SAt(S) is to find the maximum number of simultaneously satisfiable 
clauses in E(F). 

Let ni be the total number of variables used in Fi (i.e. \X l \ + \Z % \) and let uii 
be the total number of clauses in Fi. The size of F, denoted by size(F), is equal to 
Ei<i<ti(™i n i)' Given a formula E(F) specified by an L-specification F, BG{E{F)) 
denotes the bipartite graph associated with E(F). We use H [BG(E(F))\ to denote the 
L-specification of BG(E(F)). It is easy to define level-restricted l-SAt(S) formulas 



along the lines of Definition 4.4. Hence we omit this definition here. 

Example 4: Let F = (Fi (xi, X2), ^2(^3, £4), F3) be an instance of L-3SAT where each 

Fi is defined as follows: 

Fi(xi,x 2 ) = {xi +x 2 + zi) A (z 2 + z 3 ) 



F 2 {x3, x A ) = F 1 (x 3l z 4 ) A F 1 (z 4l z 5 ) A (z 4 + z 5 + x A ) 



F 3 = Fi(z 7 ,z 6 ) AF 2 (z 8 ,z 7 ) 

The formula E(F) denoted by F is (z 7 + z 6 + z J ) A {z\ + z\ ) A (z$ + z 4 + z \ ) A (z| + zf ) A 
(z 4 + z 5 + zl) A (zl + zf ) A (z A + z 5 + z 7 ). 

We now extend the definition of PL-3SAT given in |39) to define the L-PL-3SAT. 

Definition 6.3. The problem L-PL-3SAT is to decide whether the planar 3CNF 
formula E(F) specified by an L-specification F is satisfiable. The corresponding op- 
timization problem denoted by L-PL-MAX-3SAT is to find the maximum number of 
simultaneously satisfiable clauses in E(F). 

Extensions of the above definition to 1-L-PL-3SAT, 1-l-pl-max-3sat,l-pl-SAt(S) , 
l-pl-max-sat(S), I-l-pl-sat(S), I-l-pl-max-sat(S), I-fpn-pl-sat(S) 
and I-fpn-pl-max-SAt(S) are straightforward and are omitted. 

Finally we state the following PSPACE-completeness results proved in a sequel 
paper [Q . The definitions of the problems mentioned in the following theorems can 
be found in 

Theorem 6.4. The following problems are PSPACE-complete for 1 -level-restricted 
L-specified planar instances: independent set, vertex cover, partition into 
triangles and SAt(S) such that Rep(S) is the set of all finite arity Boolean rela- 
tions. 
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Theorem 6.5. The following problems are PSPACE-complete for 1-FPN- specified 

planar instances: INDEPENDENT SET, VERTEX COVER, PARTITION INTO TRIANGLES 

and SAt(S) such that Rep(S) is the set of all finite arity Boolean relations. 



7. Approximation algorithms. The hardness results in Theorems 5A and |6.5| 
motivate the study of polynomial time approximation algorithms with good perfor- 
mance guarantees for these problems. We show that several basic combinatorial prob- 



lems (including the ones in Theorems 6.4 and p.5[ ) have approximation algorithms with 
performance guarantees asymptotically equal to the best known performance guar- 
antees, when instances are specified using standard specifications. As an immediate 
corollary, most of the problems shown to have polynomial time approximation schemes 
(PTASs) in H ^l| when instances are represented using standard specifications, have 
PTASs when instances are specified either by fc-level-restricted L-specifications or 
1-FPN-specifications. 

7.1. The basic technique: Partial expansion. We outline the basic tech- 
nique behind the approximation algorithms for the 1-level-restricted L-specified prob- 
lems. Consider one of the maximization problems II in this paper. Let A be an 
approximation algorithm with performance guarantee FBEST, for II when speci- 
fied using standard specifications. Also, let T(N) denote an increasing function that 
is an upper bound on the running time of A used to solve II specified using stan- 
dard specifications of size O(N). Then, given a fixed I > 1, our approximation 
algorithm for 1-L-II takes time 0(N ■ T(N l+1 )) and has a performance guarantee of 
(!±!) . FBEST. Informally, the algorithm consists of (/ + 1) iterations. During an 
iteration i we delete^] all the explicit vertices which belong to nonterminals defined at 
level j, j = i mod (I + 1). This breaks up the given hierarchy tree into a collection 
of disjoint trees. The algorithm finds a near-optimal solution for the vertex induced 
subgraph^ defined by each small tree and outputs the union of all these solutions as 
the solution for the problem II. It is important to observe that the hierarchy tree 
can have an exponential number of nodes. Hence the deletion of nonterminals and 
the determination of near-optimal solutions for each subtree has to be done in such 
a manner so that the whole process takes only polynomial time. This is achieved by 
observing that the subtrees can be divided into n distinct equivalence classes and that 
the number of subtrees in each equivalence class can be counted in polynomial time 
in the size of the specification. 

We remark that our idea of dividing the graph into vertex (edge) disjoint sub- 
graphs is similar to the technique used by Baker || for obtaining approximation 
schemes for planar graph problems. 

7.2. MAXIMUM independent set problem for 1-level-restricted L-specified 
planar graphs. We illustrate the technique by giving a polynomial time approxi- 
mation scheme for the maximum independent set problem for 1-level-restricted 
L-specified planar graphs. The independent set problem is defined as follows. 
Given a graph G = (V, E) and a positive integer K < \V\, is there an independent 
of size K or more for 67, i.e., a subset V C V with \V'\ > K such that for each 
u, v e V (u, v) £ E ? The optimization problem called the MAXIMUM independent 
set problem (mis) requires one to find an independent set of maximum size. In 
Jicfl , we showed that given an L-specification that has edges between pins in the same 



2 For minimization problem instead of deleting the vertices in the level, we consider the vertices 
as a part of both the subtrees. 

3 For a fixed I, the size of each subgraph is polynomial in the size of the specification. 
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cell, there is a polynomial time algorithm to construct a new L-specification such that 
there is no edge between pins in the same cell. Consequently, we assume without loss 
of generality that in the given L-specification there is no edge between two pins in the 
same nonterminal. 

In the following description, we use HIS(Gi) to denote the approximate indepen- 
dent set for the graph E(Ti) obtained by our algorithm H-MIS. We also use F-MIS 
to denote the algorithm of Baker Q| for finding an approximate independent set in a 
planar graph specified using a standard specification. Before we discuss the details of 
the heuristic we define the concept of partial expansion of an L-specification. Recall 
that, for each nonterminal Gi there is a unique hierarchy tree HT(Gi) rooted at Gi. 

Definition 7.1. Let T = (Gi, G„) be an L-specification of a graph E(T). The 
partial expansion PE{G 3 i ), of the nonterminal Gi is constructed as follows: 
j = 0: PE(G J i ) = Gi— {all the explicit vertices defined in Gi} (Thus the definition 
of PE(G?) now consists of a collection of the nonterminals and pins called in the 
definition of Gi ). 

j > 1 : Repeat the following step for each nonterminal G r called by Git Insert a copy 
of PE(G^ 1 ) by identifying the I th pin of P-E(G£ _1 ) with the node in Gi that is labeled 
(v,l). (Observe that the definition of PE{G :i i ) consists of (i) explicit vertices defined 
in all the nonterminals at depth r, < r < j — 1 in HT{Gi) and (ii) a multiset of 
nonterminals Gk, such that the nonterminal G& occurs at depth j + 1 in the hierarchy 
tree HT(G l ).) 

Let Ex{PE(G 3 i )) denote the subgraph induced by the set of explicit vertices in 
the definition of PE(Gj). Also let V(E(Ti)) denote the set of vertices in E(Ti). 

Heuristic H-MIS 

• Input: A 1-level- restricted L-specification T = (Gi, G n ) of a planar graph 
G and an integer I > 1. 

• Output: An L-specification of an independent set for E(T) whose size is at 
least (j^Tj-) 2 times the size of an optimal independent set in E(T). 

• 1. For each 1 < i < I, find a near-optimal independent set in E(Ti) using 

F-MIS. 
2. For each I + 1 < i < n - 1 

(a) Compute the partial expansion PE(G\) of Gi. 

(b) Find an independent set in the subgraph Ex(PE(G\)) using heuris- 
tic F-MIS. Denote this by A\. 

(c) Let Gi i; - • • Gi p denote the multiset of nonterminals in PE(G\). 
Then the independent set for the whole graph for the iteration i 
denoted by HIS(Gi) is given by 

HIS{Gi) =A\U [j HIS{G lr ). 

l<r<p 

Remark: The explicit vertices in PE(G l i) do not have an edge to 
any of the nonterminals G^ , • • • G; . From this observation and the 
definition of hierarchical specification the independent set HIS(Gi) 
can now be calculated as follows. 

(d) 

\HIS(G t )\ = \A<\+ £ \HIS(G lr )\ 

l<r<p 
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3. For each < i < I 

(a) Compute the partial expansion PE(G] l ) of G n . 

(b) Find a near-optimal independent set of all the explicit vertices in 
PE(G l n ) using F-MIS. Denote this by A\. 

(c) Let G ni , - ■ - G np denote the multiset of nonterminals in PE{G l n ). 
The independent set for the whole graph for the iteration i, denoted 
by HISi(G n ), is given by 

HIS i (G n )=A i n U |J HIS(G nr ). 

l<r<p 

Remark: By a remark similar to one in Step 2(c) of the algorithm, 
we have the following. 

(d) 

\HISi{G n )\ = \A^\+ J2 \HIS(G nr )\. 

l<r<p 

4. The independent set HIS(G) is the largest among all the independent 
sets HISi(G n ) computed in Step 3(c). 

5. \HIS(G)\ = max \HISi(G n )\ 

0<i<l 

7.3. Analysis and Performance Guarantee. The correctness of H- MIS and 
the proof of its performance guarantee is based on the following intermediate results. 

Lemma 7.2. The set HIS(G) computed by the algorithm H-MIS in Step 4 is an 
independent set. 

Proof. We first prove that the set for 1 < i < n — 1, HIS(Gi), is an independent 
set. The proof is by induction on the depth of the hierarchy tree HT(T). 
Basis: If the depth is < I, the proof follows by the correctness of algorithm F-MIS. 
Induction: Assume that the lemma holds for all hierarchy trees of depth at most 
m> I. Consider a hierarchy tree of depth m + 1. Step 2(c) of the algorithm, computes 
a partial expansion PE(G\). This implies that the explicit vertices in PE(G\) do not 
have edges incident on the nonterminals in PE(G\). Thus, by the definition of 1- 
level-restricted L-specifications and partial expansion, it follows that the independent 
sets A\, and the sets HIS(G ir ), 1 < r < p computed in Steps 2(b) and 2(c) are 
disjoint. Also, the nonterminals in PE(G\) are at level I + 1 in HT(Gi) 1 and have an 
associated hierarchy tree of depth < m. Thus by induction hypothesis and the above 
stated observations, it follows that HIS(G m ) computed in Step 2(c) is an independent 
set. This completes the proof that 1 < i < n — 1, HIS(Gi) is an independent set. 

A similar inductive argument proves that the set HISi{G n ) computed in each 
iteration of Step 3(c) is also an independent set. By Step 4, we have that HIS(G) is 
an independent set. □ 

Lemma 7.3. 

1. In each iteration i, I + 1 < i < n — 1, of Step 2 of algorithm H-MIS, all 
the explicit vertices in nonterminals at levels j = I mod (/ + 1) in the hierarchy tree 
HT(Gi) are deleted. 

2. In each iteration i of Step 3 of algorithm H-MIS ; all the explicit vertices in 
nonterminals at levels j — i mod (/ + 1) in the hierarchy tree HT(G n ) are deleted. 

Proof. 

Proof of Part 1: Induction on the depth of the hierarchy tree associated with d. 
Basis: If the depth is l + l, the proof follows directly by Step 1 and the definition of 
partial expansion. 

Induction: Assume that the lemma holds for all hierarchy trees of depth at most 
m > (I + 1). Consider a hierarchy tree of depth m + 1. Step 2(c) of the algorithm, 
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computes the partial expansion PE(G\). This implies that all the explicit vertices at 
level I in the hierarchy tree HT(Gi) were deleted. Each nonterminal occurring in the 
definition of PE(G\) is at level 1+ 1 in HT(Gi), and has an associated hierarchy tree 
of depth < to. The proof now follows by induction hypothesis. 

Proof of Part 2: Consider a hierarchy tree HT(G n ). In iteration i of Step 3 we 
compute PE(G\). This removes all the explicit vertices defined in nonterminals at 
level i. Also, by the definition of partial expansion it follows that all explicit vertices 
defined in nonterminals at levels 1 to i appear explicitly in the partially expanded 
graph. Therefore, the partially expanded graph now has nonterminals defined at level 
i + 1 in the hierarchy tree HT(G n ). The theorem now follows as a consequence of 
Part 1 of the theorem. □ 

Given the decomposition of E(T) into a forest (as a result of removing explicit 
vertices, in nonterminals at levels j — i mod (I + 1) in the hierarchy tree HT(G n )) 
we can associate a hierarchy tree with each of the subgraphs in the forest. Each such 
tree is a subtree of the original hierarchy tree HT(T). Label each subtree by the type 
of nonterminal that is the root of the subtree. The proof of the following lemma is 
straightforward . 

Lemma 7.4. 

1. During each iteration i of Step 3 of the algorithm H-MIS, the root of each 
subtree is labeled by one of the elements of the set {G\, ■ • • , G n _i}. 

2. For 1 < i < n, let H\, . . . , H^. be the set of graphs corresponding to the 
subtrees labeled Gj. Then for each i the graphs H\, . . . , H*. are isomorphic. 

7.3.1. At Least One Good Iteration Exists. Next we prove that, at least 
one iteration of Step 3 has the property that the number of nodes of an optimal 
independent set that are deleted is a small fraction of the optimal independent set. 
Let Fi denote the set of vertices obtained by d eleting the explicit nodes in iteration 



i in Step 3 of algorithm H-MIS. By Lemma 7.3 it follows that for each iteration i we 
did not consider the explicit vertices in levels , ji 2 • ■ ■ ji such that 1 < i p < n and 
ji = i mod (£ + 1), 1 < q < p. Let Si, < i < I, be the set of vertices not considered 
in iteration i of Step 3. Let IS{G n ) denote an optimum independent set in the graph 
E(T). Let IS pt(Si) denote the nodes in Si included in the maximum independent 
set IS(G n ). 
Lemma 7.5. 

max \IS(Fi)\ > -J—\IS(G n )\ 
0<t<l [i -\- 1) 



Proof. By Lemma 7.3 and the algorithm H-MIS, it follows that 
S, n Sj = 0, St = V{E(T)), and 

\IS opt (S )\ + llS^tiSx)] + ■■■ + \IS opt (Si)\ - \IS(G n )\ 

Therefore, 

min |JS op tOSj)| < \IS(G n )\/(l + l) 

0<i<l 



max 1/5(^)1 > |/5(G n )| - min \IS opt (S t )\ > -—\IS(G n )\. 

0<i<l 0<i<( u + 1) 

□ 



16 M.V. MARATHE, H.B. HUNT, III, R.E. STEARNS AND V. RADHAKRISHNAN 



7.3.2. Performance Guarantee and running time. We now prove that the 
above algorithm computes a near-optimal independent set. Given any e > 0, for some 
choice of positive integer I such that (j^j-) 2 > (1 — e), we show that algorithm H-MIS 
computes an independent set whose size is at least (1 — e) times the size of an optimal 
independent set. We first recall a similar lemma in []3| for planar graphs specified 
using standard specifications. 

Theorem 7.6. J|/ For all fixed 1>1, given a planar graph G there is linear time 
algorithm that computes an independent set FIS(G) such that \FIS(G)\ > (jtt[) • 
\IS(G)\, where IS(G) denotes a maximum independent set in G. 

Lemma 7.7. \HIS l (G n )\ > (j^) ■ \IS(F t )\. 

Proof. Induction on the number of nonterminals in the definition of T. The base 
case is fairly straightforward. Consider the induction step. By the definition of partial 
expansion it follows that, 

\IS(F l )\ = \IS(Ex(PE(Gl l ))\+ J2 \IS(PE(G nr ))\. 

l<r<p 

From Step 3(c) of the algorithm H-MIS we also know that 



\HISi(G n )\ = K| + \HIS{G nr )\. 

l<r<p 

From the induction hypothesis and Theorem |7.6| it follows that 

Kl > {j^)-\IS(Ex{PE{Gi))\ and 

\HIS(G nr )\ > (j^)-\IS(PE(G nr ))\. 
The lemma now follows. □ 

Theorem 7.8. \HIS(G)\ > (jtt) 2 • \IS(G)\. 



Proof. Follows from Lemma 7.5 and repeated application of Lemma 7.7. □ 
Theorem 7.9. Let T be an L- specification with vertex number N. Given any 
e > 0, let I > 1 be an integer such that (-fp[) 2 > (1 — e)- Then the approximation 
algorithm H-MIS runs in time 0(N l+2 ) and finds an independent set in E(T) that 
is at least (]zp[) 2 times the size of an optimal independent set in E(T). 



Proof. The performance guarantee follows by Theorem 7.8. Therefore we only 
prove the claimed time bounds. 

First consider Step 1. Note that by Euler's formula, the number of edges in a 
planar graph with 0(N l ) vertices is also 0(N l ). Thus, the size of the graphs E(Ti), 
1 < i < I is 0(N l ). Hence the time required to compute the partial expansion is 
0(N l ). By Theorem |7~6| , the time needed to compute an independent set in E(Ti) is 
0(N l ). Thus the total running time of Step 1 is 0(N l ). 

Next consider each iteration of Step 2 of the algorithm H-MIS. Step 2(a) takes 
time 0(N l+1 ) since the size of the graph PE(G\) can be 0{N l+1 ). By Theorem 



7.6, the time needed for executing Step 2(b) is 0(N ), since the number of nodes 



in Ex(PE(G\)) can be 0(N l ). By Lemma |7J, Step 2(c) and 2(d) together take 
time O(N). Therefore the total running time for executing one iteration of Step 2 is 
0{N l+1 ). Thus the total running time of Step 2 is nO{N l+1 ) = 0{N l + 2 ). 

A similar calculation shows that the total time needed to execute one iteration of 
Step 3 is 0(N l+1 ). Thus the total time needed to execute Step 3 is (I + l)0(N l+1 ) = 
0{N l+1 ). 

Thus the total running time of the algorithm is 0(N l+2 ). □ 

7.4. L-Specification of the solution and the query problem. In §^^, we 
showed how to solve the size problem for 1-l-mis. We now discuss the construction 
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problem. As noted in §2.1 our algorithms for the four variants of the problem apply 
to the same independent set HIS(G). 

The L-specification of the solution can be easily constructed by slightly modify- 
ing the algorithm H-MIS as follows. Consider the iteration i of Step 3 which gives 
the maximum independent set. Denote the iteration by i*. The L-specification H 
of the solution consists of nonterminals Hi, • • • , H n . For 1 < j ' < n the explicit ver- 
tices of H,j are the explicit vertices in PE{G^) that are in the independent set. If 

PE(G l j) calls nonterminals Gj 1 , ■ ■ ■ ,Gj m then the nonterminal Hj calls the nonter- 
minals Hj 1 , ■ ■ • , Hj m . Observe that some of the nonterminals Hi may be redundant 
and these can removed from the final specification. Given the L-specification of the 
solution, the query problem be easily solved by examining if the given vertex oc- 
curs in the set of nodes specified by the L-specification of the solution. Given an 
L-specification of the solution, we can solve the output problem as follows: We tra- 
verse the hierarchy tree associated with H in a depth first manner and output the 
vertices in the nonterminals visited during the traversal. 

Observe that the only place we used planarity was to obtain a near-optimal so- 
lution for the maximum independent set problem for each partially expanded graph. 



In §7.7 we use this observation to compute near-optimal solutions for problems for 



arbitrary 1-level-restricted L-specified graphs. 

7.5. Other L-specified planar problems . Our technique can be applied to 
obtain efficient approximation algorithms for the following additional optimization 
problems: minimum vertex cover, maximum partition into triangles, mini- 
mum edge dominating set, MAXIMUM CUT and max SAt(S) for any finite set of 
finite set of finite arity Boolean relations S. The basic idea behind devising approxi- 
mation schemes for these problems is similar to the ideas used to solve the maximum 
independent set problem. Therefore, we only briefly discuss the method for mini- 
mum VERTEX COVER and MAX SAT(S). 

(1) MINIMUM vertex COVER:. Given a graph G = (V,E) and a positive integer 
K < \V\, is there a vertex cover of size K or less for G, i.e., a subset V C V with 
\V'\ < K such that for each edge (u,v) e E either u or v belongs to V' ? The 
optimization problem requires one to find a vertex cover of minimum size. 

In order to approximate the 1-l-pl-MINIMUM vertex cover problem we do the 
following. Given an e, we choose an / such that (^f-) 2 < (1 + e). Next, we modify the 
definition of partial expansion so that instead of deleting the explicit vertices at levels 
(/ + 1) apart, we consider them in both sides of the partition. For each < i < I, 
the algorithm finds a near-optimal solution for the overlapping planar graphs induced 
by explicit vertices in levels (jl + i) to ((j + 1)1 + i), for j > 0. The algorithm 
picks the best among all the vertex covers obtained for the different values of i. Let 
OPT(G) denote an optimal vertex cover for G. The following lemma points out that 
the solution obtained is at most (^f-) 2 times the optimal vertex cover. The proof of 
the lemma follows the same general argument given for the MAXIMUM independent 
SET problem. 

Lemma 7.10. The size of the vertex cover obtained is no more than 

( l -±lf\OPT(G)\ 



Proof. Consider an optimal solution OPT(G) to the vertex cover problem. Then 
for some < t < I, at most \OPT(G)\/l nodes in OPT{G) are in levels congruent to 
t mod (I). Consider the iteration when the planar graphs are obtained by overlapping 
at levels congruent to t mod (I). Hence the size of an optimal vertex cover in this iter- 
ation is (\OPT(G) \ + \OPT(G)\/V). Now applying the known approximation scheme 
H for computing a near-optimal vertex cover for each of smaller subgraphs, we obtain 



18 M.V. MARATHE, H.B. HUNT, III, R.E. STEARNS AND V. RADHAKRISHNAN 



Variables in level j included in this 
subgraph 





Clauses removed 
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Variables in level (j+1) 
included in this subgraph. 



Fig. 7.1. Basic idea behind the approximation algorithm for I-l-max-pl-sat(S) . The black dots 
represent variables and the ellipses denote clauses. The figure depicts the set of clauses to be deleted 
and the redistribution of the variables. 



a near-optimal vertex cover for the whole graph for iteration t. The size of the vertex 
cover obtained in this iteration is no more than (\OPT(G)\ + \0 PT '(G) | / 'Z)^±I . The 
reason is that the explicit vertices in the overlapping levels are counted twice and 
the near-optimal vertex cover heuristic yields a vertex cover of size (Z + 1)/Z times 
the optimal vertex cover for each subgraph. Since the heuristic picks the minimum 
vertex over all values of i, it follows that the size of the vertex cover produced by the 
heuristic is no more than ( l -±±) 2 \OPT{G)\. □ 

(2) max SAt(S):. In the following, we will assume that an instance F of 1-L- 
pl-max-SAt(S) is specified by H[BG{E{F))] (i.e the specification of the associated 
bipartite graph.). The basic idea behind the approximation schemes for 1-l-max- 
pl-SAt(S) is as follows: For each i, < i < 21 in increments of 2, we remove the 
explicitly defined clauses which are in levels j and j + 1, such that j = i mod (Z + 1). 
This breaks the bipartite graph into a number of smaller bipartite graphs such that 
the formulas they denote do not share any variables or clauses. It is not difficult 
to modify the defin ition of partial expansion to obtain a decomposition as described 



above. Figure 7.1 shows how the variables in levels j and j + 1 are redistributed. 
As in the case of MAXIMUM independent set problem, it is easy to see that there 
exists an iteration t, < t < 21, such that at most jjjpfj clauses in OPT are deleted. 

Next, by the results in pi| the problem can be solved near-optimally for each smaller 
subformulas. The union of the clauses satisfied for each small formula constitutes a 
solution for a given value of i. We pick the best solution for different values of i. This 
ensures that the best assignment to the variables over all values of i is at least (j4j) 2 
of an optimal assignment to the variables of the I-L-PL-MAX-SAt(S) instance. 

7.6. Extension to fc-level-restricted instances. The technique used to solve 
various problems for 1-level-restricted L-specifications can be generalized to solve 
problems specified using fc-level-restricted L-specifications. We only point out the 
essential differences. Again, for the purposes of illustration consider the problem 
fc-L-PL-MiS. First note that we need to extend the definition of partial expansion 
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so that we delete the explicit vertices in nonterminals at k consecutive levels. This 
implies that the time to compute PE{G\), 1 < i < n - 1 is 0(N l+k ). The rest of the 
algorithm follows the same outline as that of H-MIS. The pro of of correctness and 
the performance guarantee also follow similar arguments as in § |7.3| . Thus the total 
running time of the algorithm is 0(N k+l+1 ) and its performance guarantee is (^f-) 2 - 
Hence we have the following theorem. 

Theorem 7.11. For any fixed k > 1, there are polynomial time approximation 
schemes for the problems MAXIMUM INDEPENDENT SET, MINIMUM VERTEX COVER, 
MINIMUM EDGE DOMINATING SET, MAXIMUM PARTITION INTO TRIANGLES AND MAX- 
IMUM CUT, and max SAt(S), for each finite set of finite arity Boolean relations S. 
when restricted to planar instances specified using k-level-restricted L- specifications. 

7.7. Ext ension to level restricted arbitrary instances. Our results in §[7 



through §7.6 can be extended for problems on arbitrary graphs specified using k- 
level-restr ictc d L-specifications. To do this, observe that to obtain the results in § |7.2| 
through § |7.6| we used planarity only to obtain approximation schemes for smaller 
subgraphs (formulas) obtained as a result of partial expansion. If the graphs were 
not planar we could use the best known approximation algorithms for solving the 
problem near-optimally and in turn get a performance guarantee which reflects this 
bound. For example, consider the problem 1-l-max-2sat. Let e > be the required 
performance guarantee. I > 1 is an integer satisfying the inequality > (1 — 

e). For the problem MAX-2SAT, the recent work of Goemans and Williamson JlJ] 
provides an approximation algorithm with performance guarantee of 1.137. Using 
their algorithm as a subroutine to solve the small MAX-2SAT instances obtained as 
a result of partial expansion, we can devise an approximation algorithm for 1-L- 
MAX-2SAT with performance guarantee (^-) 1.137. A similar idea applies to other 
optimization problems considered. Again, it is easy to generalize our results for k- 
level-restricted L-specifications. Thus we have the following theorem. 

Let II be one of the problems: maximum independent set, MINIMUM vertex 

COVER, MINIMUM EDGE DOMINATING SET, MAXIMUM PARTITION INTO TRIANGLES, 

maximum CUT and MAX-SAt(S), for finite set of Boolean relations S, such that Rep(S) 
is the set of all finite arity Boolean relations^]. 

Theorem 7.12. For all fixed k > 1, e > and for all of the problems LI, 
there are polynomial time approximation algorithms with performance guarantee^ (1 + 
e) • FBESTu for problems II, when specified using k-level-restricted L-specifications. 
Here FBESTu denotes the best known performance guarantee of an algorithm for the 
problem LI for instances specified using standard specifications. 

Using the results of Arora et al.j|, Bellare et. al. Q and our results in |2^] we 
get the following theorem. 

Theorem 7.13. Unless P = NP, the problems II, when specified using k-level- 
restricted L-specifications, do not have polynomial time approximation schemes. 

7.8. Approximation algorithms for 1-FPN-sp ecifi ed prob lem s. Next, we 
briefly discuss how to extend our ideas developed in § |7.2| through §|7.7| in order to 
devise approximation schemes for several PSPACE-hard problems for 1-FPN-specified 
instances. 

The basic idea is simple. Once again we illustrate our ideas by describing our 
approximation algorithm for the problem 1-fpn-pl-mis. Given a 1-FPN-specification 
L = (G(V, E),m) of a planar graph G m and an e > 0, we find the corresponding 
integer I that satisfies the inequality (j^i) 2 > (1 — e). For < i < I, we remove the 
vertices placed at the lattice points j such that j • = i mod(l + 1). This partitions the 
graph G m into a number of smaller disjoint subgraphs, each induced by I consecutive 
lattice points. 



4 Actually our easiness results hold for all finite set of finite arity Boolean relations S. 

5 For the sake of uniformity we assume that the performance guarantee is > 1. 
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Fig. 7.2. A schematic diagram showing the vertices to be removed in each iteration i while 
computing a near-optimal independent set for 1-FPN-specified planar graphs. In our example i = 3, 
I + 1 = 4, and m = 9. Each box represents a copy of the vertices in the original static graph. The 
shaded area represents the vertices that are removed. 



Specifically, for a given i, let l l p — max{0, (p — 1)(Z + 1) + (i + 1)} and r p = 
mm{m,p(l + 1) + (i — 1)}, where < p < fcj. Here i, = |" m n+i) 1- ^ e su ^~ 
graph induced by vertices v(j p ), where l p < j p < r*, be denoted by H(JL,r p ). For 
a given e > 0, the graphs H(l p , r p ) are linear in the size of T. Figure [7^ shows a 
schematic diagram of the vertices removed in a given iteration i. Next, we solve the 
mis problem near-optimally on each of th e su bgraphs. This can be done by using the 



linear time algorithm stated in Theorem 7.6. The union of these independent sets is 



the independent set obtained in iteration i. The heuristic simply picks up the largest 
independent set obtained over alH + 1 iterations. By arg um ent s sim ilar to the ones 



we presented for approximating 1-l-pl-mis (Subsections [7^ to 7.4), it follows that 
the approximation algorithm has a performance guarantee of (^-) 2 • 

We note the following important point. If a near-optimal independent set were 
to be obtained for each subgraph H (l^, r 1 ^) , we would take an exponential amount 
of time in each iteration i. This is because p = 0(m). Hence we can not afford to 
solve the problem explicitly for each subgraph. But observe that each iteration i the 
subgraphs H(l p ,r p ), 1 < p < (7+1)^ 1 ~ 1 are isomorphic. Hence we need to solve the 
MIS problem for the graphs rj), H(l\, r\) and H(l\.,r\.), where U = \ ™^£=*L \ 
Let IS(H(l p , r l p )) denote the independent set obtained by the heuristic for the graph 
H(l pl r p )). Furthermore, let the approximate maximum independent set for the whole 
graph for a given iteration i be denoted by IS(G m (i)). Then the size of IS(G m (i)) 
is given by the following equation: 

\IS(G m (i))\ = \IS(H(liri))\ + L m ( ~ |7) ^ }\ IS ( H ( l l4))\ + \IS(H(li,ri))\ 

This completes the discussion of the approximatio n alg orithm for 1 -fpn-pl-mis. 
By combining the above arguments along with those in § [7.2| through §7.7, we can show 
that several other optimization problems can be approximated in a similar fashion. 
Again it is easy to see that the technique extends to problems for arbitrary instances 
and also to problems for instances specified using fc-narrow 1-FPN-specifications. 
Thus we have the following theorem. 

Theorem 7.14. For all fixed k > 1, e > and for all of the problems II stated in 
§]7. 9| there are polynomial time approximation algorithms with performance guarantee^ 
(1 + e) • FBESTn for problems Ti, when specified using k-level-restricted 1-FPN- 
specifications. Here FBESTu denotes the best known performance guarantee of an 
algorithm for the problem H for instances specified using standard specifications. 

Observe that the technique used to devise approximation algorithms for problems 
restricted to fc-narrow 1-FPN-specified instances is very similar to the technique used 
to devise approximation algorithms for /c-level-restricted L-specified problems. But 
there are two important differences in the details of the algorithms. 



For the sake of uniformity we assume that the performance guarantee is > 1. 
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1. In case of algorithms for L-specified problems, the number of equivalence 
classes is 0(n) where n is the number of nonterminals. In contrast, the number of 
equivalence classes in case of algorithms for 1-FPN-specified problems is only O(l). 

2. The size of the subgraphs for which the problem is solved near-optimally 
also differs significantly. Specifically, the number of explicit vertices in PE(G\) can 
be 0(N l ). Moreover the time required to compute PE(G\) can be 0(N l+k ). In 
contrast, the number of explicit vertices in each H(l p , r p ) is only O(N) and the time 
required to construct each H(l l p , r p ) is only O(N). In both cases we use N to be the 
vertex number of the respective specifications T ( N can be 0(size(T)). 

These important differences allow us to devise linear time approximation schemes 
for 1-FPN-specified problems. 

8. Conclusions. 

8.1. Summary. We have investigated the polynomial time approximability of 
several PSPACE-hard optimization problems for both L- and 1-FPN-specified in- 
stances. A general approach was given to obtain polynomial time approximation 
schemes for several PSPACE-hard optimization problems for planar graphs specified 
using fc-level-restricted L- or 1-FPN-specifications. We believe that the partial ex- 
pansion technique can be used to obtain efficient approximations for other problems 
specified using L- or 1-FPN-specifications as well as for problems specified using other 
succinct specifications. 

In an accompanying paper fii| , we investigate the decision complexity of various 
combinatorial problems specified using various kinds of L-specifications and 1-FPN- 
specifications. There we give a general method to obtain PSPACE-hard lower bounds 
for such problems including the ones discussed here. 

8.2. Open Problems. We conclude with a list of open problems for future 
research. 

1. Can we use the concept of Probabilistically Checkable Debate systems ||, |t| 
to prove non-approximability results for problems specified using arbitrary (not level- 
restricted) L-specifications ? 

Recently, Agarwal and Condon OJ have partially answered this question by showing 
that unless P = PSPACE, there is no polynomial time approximation scheme for the 
problem L-MAX-3SAT. The result was proved by using the characterization of PSPACE 



in terms of random debate systems. In 
L-MAX-SAt(S) such that Rep(S) denotesT 



we extended their result to hold for any 
re set of all finite arity Boolean relations. 



2. Recently, several researchers have considered logical definability of a number 
of optimization problems and defined appropriate classes such as MAX SNP MAX IT 
MAX NP and MAX #P (cf. pL M ~ 



J ) . All these researchers have assumed the 
the input is specified using standarcTspecifications. What happens if the instances 
(finite or infinite) are specified succinctly ? 

Some work has been done along these lines by Hirst and Harel |l7j . Specifically, they 
considered infinite recursive versions of several NP optimization problems. They prove 
that some problems become highly undecidable (in terms of Turing degrees) while 
others remain on low levels of arithmetic hierarchy. As a corollary of their results 
they provide a method for proving (finitary) problems to be outside the syntactic 
class MAX NP and hence outside MAX SNP. 
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